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ABSTRACT 


Recent observations of cooling flare loops indicate that the differ- 
ential emission measure, Q, of flare coronal and transition region plasma, 

T > 10® K, has a much steeper dependence on temperature than in non-flare 
regions. This result is not compatible with models in which conduction 
dominates the cooling; hence, we investigate maJels in which radiation 
dominates. We find that the radiative models predict Q « T^"^’ , where Jl 
measures the dependence of the radiative loss coefficient on temperature, 
A(T) ~ T~^', We conclude that the radiative models are also incapable of 
explaining the observations and suggest that large mass motions (velocities 
of the order of the sound speed) may be required. 



I. INTRODUCTION 


The high temperature coronal plasma produced by a flare is believed 
to cool primarily by radiation and conduction to the chromosphere. Since 
conductive losses increase rapidly with teraperaturcj oeveral authors con- 
cluded that conduction dominates radiation, at least during the initial 
phase when T > 10^ K, e.g. , Culhane et al. (1970), Moore and Datlowe (1975). 
However, radiative losses are a strong function of density, and recent EUV 
and soft x-ray observations have tended to yield increasingly higher density 
estimates; n > lO'^ cm"^ at T s lO’ K, (see Moore et al. 1979). Hence, it 
appears likely that radiative cooling may dominate in some flares even 
during Che initial cooling phase. 

Strong evidence against Che dominance of conduction has been obtained 
from observations of relative line intensities of coronal and transition 
region lines, i.e. 10® K < T < 10^ K, Given a model for the cooling, the 
temperature and density profiles and, hence, the emission measure pro- 
file (differential emission measure) can be determined. This may then be 
used to predict relative line intensities which can be compared with the 
observations. Conductive dominated models have been investigated in detail 
by Antiochos and Sturrock (1976, 1978). For a plane-parallel isobaric 
atmosphere the differential emission measure, 


q(T) 


= A n'^ 


1 ^ 
T 3S 




(i) 


is related directly to the conductive heat flux. 
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F 


c 



( 2 ) 


where; A Is the area of the emitting region; n is the electron density; 
9T 

is the temperature gradient in the vertical direction; and the con- 
ductivity K is given by Spitzer (1962) : 


K = 10 -^ T®/' 


(3) 


Combining (1) - (3) yields: 


10-®Ap^ 

Q(T) ["^ (A) 


where the equation of state, 


p = 2ktiT 


(5) 


has been used. 

Antiochos and Sturrock find that for static conductive cooling, « 
const whereas for evaporative conduction “ T. Defining 5 as the steep- 
ness of the differential emission measure, Q(T) « equation (4) implies 
that 6 = 1.5 for static conduction and 6 » ,5 for evaporative. Note that 
these results apply only to a single coronal loop (of constant cross-section), 
but that the observed emission measure may be due to the contributions of 
several loops. However, if all the loops are conductively cooling then 
the observed value of 6 must be < 1.5 since no individual loop can have 


6 > 1.5. 
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The conductive models are conoistant with EUV observations of quiet 
and active regions which Indicate that 6 < 1,5, e.g. , Jordan (1976), Withbroe 
(1977), However, the models disagree with recent flare observations which 
imply a very steep dependence of emission measure on temperature. Dere 
et al, (1977), Dare and Cook (1979), and Widing (1979) find that in three 
flares 6 > 3,0 for plasma at 10®*^ < T 5 lo’ . Underwood et al. (1978) carried 
out a detailed comparison of observed line intensities for the 1973 August 9 
flare with predictions by the conductive models. By comparing intensities 
directly they avoided any difficulties with deriving the differential 
emission measure from the data (Craig and Brown 1976, Underwood and McKenzie 
1977). They also found that only models with a large 6, > 3.0, could be 
compatible with the data, and that the conductive models were definitely 
in disagreement. It is difficult to reconcile a large 6 with conduct:tve 
cooling since equation (4) implies that fur 6 > 1,5 the heat flux actually 
increases as it passes through the transition region. 

Dere and Cook (1979) have attempted to explain their observations 
as due to the effect of a collection of loops. However, as pointed out 
previously, no sum of loops can have a differential emission measure 
steeper than, say 6 =1.5, unless at least one loop has this value or 
greater. Dere and Cook assume an infinite steepness, i.e, , they assume 
that each loop is perfertly isothermal and, hence, the differential emission 
measure of each loop :,s a delta function; but, such a procedure is physically 
meaningless. Clearly, the temperatui-e profile in each loop must be a 
continuous function from its maximum value in the corona to its value in 
the chromosphere, ~ lO** K, and must be determined by the physical processes 
in the plasma, e.g,, conduction and radiation. The point of the observations 
in that they impose a restriction on the possible form of the temperature 
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profile and, hence, on the physics of the loop. Dere and Cook have failed 
to understand this point and have ignored the whole issue by simply assuming 
an impossible temperature profile. 

Underwood et. al. proposed that a possible explanation for the steep 
profile is that it is due to radiative dominated cooling. The radiative 
losses are known to be a strong function of temperature and, hence, may 
be expected to result in the preferential depletion of low temperature 
plasma over that at higher temperatures, In addition, these authors observe 
large downward velocities at temperatures where che radiative losses peak, 

T ~ 10® K. Downward velocities are also difficult to reconcile with con- 
ductive dominated cooling, but they are compatible with the idea that thermal 
instability creates large pressure gradients which subsequently generate 
mass motions, (Antiochos 1979). 

From the discussion above it appears likely chat conductive models 
are invalid, at least, for some flar_‘S. Therefore, in the next section 
we develop models for the cooling of coronal flare plasma in which radiation 
dominates and, in particular, we investigate the emission measure profiles 


that such models predict. 



II. MODELS 


2.1 Approximations 

Assuming a plane-parallel geometry (a loop of constant crosa-sectlon) , 
and neglecting gravity, the heat equation can be expressed as: 

( 3 / 2 P + i/zpv^) + ^ (i/zpv^ + s/zpv + Fg) = -n^A(T) (6) 


where: A(T) is the radiative loss coefficient, e.g., Cox and Tucker 1969; 

p is the mass density; and it is assumed that all flare heating has ended 
so that there is no energy source term present. (The steady-state coronal 
heating rate is negligible compared to the energy loss rate of flare plasma.) 
Simultaneous with (6) the equations of continuity and momentum must be 
solved. However, the resulting system of equations is sufficiently complex 
that analytic general solutions are not possible. In addition, numerical 
simulation is impractical because an extremely fine spatial grid is required 
to resolve the large temperature gradients expected at lower temperature. 
Hence, we investigate only certain limiting cases to the complete set of 
equations. 

There are three physical time scales that are relevant to the model; 
the radiative cooling time 


(T) 


3/2p 

n^A(T) 


(7) 


the conductive cooling time, 


- 6 - 


and the speed of sound propagation time 


(T) H 


H(T) 

C(T) 


where H(T) is the temperature scale height, 


H(T) 


3 ^ ^ -1 

T 3a 


( 9 ) 


( 10 ) 


and C(T) is the speed of sound > We assume throughout that radiation dominates 
the cooling, 


(T) 
Tc (T) 


< < 1 


( 11 ) 


for T in the range of interest, and hence, the heat flux in equation (6) 
can be neglected. Two limiting cases can now be defined for the evolution 
of the flare plasma: 


the static 


tr (T) 

Ts (T) 


< < 1 


(12a) 


and the isobaric Zr .... ( . T . ) . > > i 

T g 


(12b) 


In the static case (12a), the radiative cooling is so rapid that even though 
large pressure gradients are created due to the different cooling rates of 
plasma at different temperatures, there is insufficient time for large 
velocities to be generated. Therefore, if the plasma is initially static 
it will remain so during the cooling phase. In. the isobaric case (12b), 
the speed of sound is so large that mass motions cancel the build up of 
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pressurc gradients so that the plasma is approximately isobaric aitrt all 
velocities remain small compared to the speed of sound. 

Conditions (11) and (12a) or (12b) apply to Clares with a high coronal 
density os in the 1973 August 9 event, A?jsumlng that fi(T) » 10“*^’® for 
T = 10’ K yields: 

rp7 / 2 

— = 10’ ® (13) 

(nH)^ 


T T** ^ * 

— » lo” ~ (14) 

(nU) 

Hence, for T = lO’ K and H « 10®'® cm, radiation dominates of n > lO’®'® 

cm”^ . The isobaric case is valid if lO’*'® > n > lO’ “ ‘ ® , whereas the static 

case Itolds of n > lo’^'® cm However, these numbers are highly variable 

since equations (13) and (14) indicate that T /T and T /t are sensitive 

r c r s 

functions of T, n and H, which are not accurately measured. In addition, 
the plasma parameters vary with time, therefore, each of the radiative 
cases may be valid only over a limited range of the evolution. 

The important points of (13) and (14) are that conditions (il) and 
(12) are mutually consistant and that both of the radiative cases are 
possible in view of the uncertainties in T, n and H. Even though neither 
of the radiative models is complete, much can be gained by examining these 
simple models before attempting extensive numerical simulations; in par- 
ticular, it will be possible to determine whether radiative cooling can 
produce a steep emission measure profile. 


2.2 Static Radiative Cooling 

Along with assumptions (11) and (12a) initial conditions must be specified 


on the plasma profiles, For simplicity, and in ordar to be consistent 
with assumptions (12a), we assume that the plasma ip initially static 
and Isobaric, 


v(s,0) = 0 and (s»0) “ 0 


( 15 ) 


Rather than specifying the Initial temperature profile, it is equivalent 
to specify the initial emission measure profile. This is assumed to be 
given by the static conductive model of Antlochos and Sturrock 1976: 


Qo<To) 


(r v) 




( 16 ) 


where V is the volume of the loop; n^ and are the density and temperature 
at the loop apex (the minimum and maximum value resp.) ; and the subscript 
•'o" indicates that the variable is to be evaluated at the time t = 0. 

By using (16) as the Initial profile, the physical situation that is being 
investigated corresponds to that of a flare loop that is initially cooling 
by conducticii, but then radiation cooling begins to dominate and, hence, 
the temperature profile changes. The evolution of the profile and of the 
differential emission measure is calculated below. 

Using (11), (12a) and (15) the plasma equations reduce to: 



-9 


and the equation of state, (5). 

Asauming a simple power law form for the radiative losses: 

A(T) « LT“^ ( 18 ) 


where L and are constants, equation (17) can be solved analytically to 
yield T(a,t): 


T(s,t) « T„(s) 


11 - 


(1 + A) t' 
T,(T^) . 


1/ft+l 


(19) 


where T^(T^) is the radiative cooling time, equation (7), for plasma at 
temperature T^t 


T (T ) 
r 0 



(T^(s))“^~^ 


( 20 ) 


The differential emission measure, equation (1), can now be determined. 
Expressed as a function of T and t it is given by; 

Qq(To) T*'+i 

Q(T,t) - — (21) 

where is also considered as a function to T and t, and is obtained 
by solving equation (19) for TQ(T,t). 

In Figure (1) , Q is plotted as a function of T/T^^ at the times t = 

0 and tj and for the value of i = 1.5. This is a somewhat larger value 
for Si than is currently accepted (e.g., Rosner et al. 1978), however, we 
selected it so that the effect of radiative cooling could be clearly 
visible in the figure. The time tj was chosen so that the maximum temper- 
ature has decreased by 10% from its initial value, i.e., T(Sj^, tj) = 
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.9 T (s ) ■• ,9 T , where s is the position of the loop npex. 

As noted previously, the initial profile has a value of 6 ■ 1.5, 

except nsar the singular point where 6 becomes infinite. This 

singularity has no physical significance (it does not mean that the lutcn;^ity 
of lines formed at this temperature ore infinite) ; it simply reflects the 
mathematical fact that at a temperature maximum the temperature gradient 
vanishes and, hence, the transformation from s to T implicit in the defi- 
nition of Q, equation (1), is no longer valid. The emission measure Q 
cannot b<, used to calculate line intensities near these temperatures, one 
must use the temperature and density profiles directly. If this is done, 
one notes no significant deviation from the ^ relation, e.g,, 

Underwood et al. (1978), ■ 

It 'R- f-.ident from the figure that radiative cooling can produce a 
steep emission measure profile, at least, for the value of I = 1.5. The 
curve for Q(T,tj) indicates that 6(tj) = 2.5 except near the singular 
point T .9 T^. This result can be obtained directly from equations 
(19) and (21). For T(s, tj ) /T^(s) « 1, i.e., after a significant amount 
of cooling has occurred, equation (19) Implies that T^(T,t ) is approximately 
independent of T and, thus, equation (21) implies that Q(T,tj) « . 

On the other hand for T(s, t, ) /T^(s) ^ 1, i.e., for material that has cooled 
very little, then equations (19) and (21) imply that Q * as expected. 
Although at t = tj the hottest plasma has a temperature only 10% less than 
its initial value, material at lower temperatures have cooled relatively 
much more due to the strong temperature dependence of the radiative cooling 
time, equation (20). For example, ^the plasma with temperature T(tj) = 

10“* T had an initial temperature of T !» .7 T . Thus, the condition 
that T(s,tj)/T^(s) « 1 is valid for all temperatures except, again. 
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for a negligibly small Interval near the maximum. We conclude, therefore, 
that for static radiative cooling the steepness of the differential emission 
measure is given by 6 = S. + 1. 


2.3 Isobaric Radiative Cooling 

Under assumptions (11) and (12b) the plasma equations become; 


^ + ^ = 0 
3t 8s 


(22a) 


in = 

8s 


0 


(22b) 


and 


1 in 

2 8t 


+ |te.VA(T) 


(22c) 


Using the approximation for the radiative losses (18) , the set above may 
be combined to yield a single equation: 


1 dp r 

5p dt . 


S 3 yi bp 

- 2(Jl + 2) y 

8s^. 


lOk^ Ly 


.8s/ 




a + 2) 


8*y 


8y 8^y 8^ y 8y 


8s^ J 8t 8s^ 8t8s 8s 


= 0 


(23) 


where y is defined as. 


y = 


(24) 


Equation (23) is too complex to solve generally for arbitrary Initial 
conditions; however, particular solutions may be obtained by separation 
of variables. Letting, 
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y(s,t) = T(0,0) 0(t) ?(8), 


(24) 


and 


P(t) = 


(25) 


equadion (23) reduces toJ 


5 T d0 / ( 


T 0 d(j) 
1 0 ^ 

u- 



1+ . 

3 (}) dt \ 

d^ da ! 

3 (|)^ dt 

d5^ 

d^e' 


d=^e 1 dc’^ 

(22 + 7) — - 2(2 + 2) 



(2 + 2) + 

ds 

ds^ „ 


ds^ ? ds 


where 

3/2 P 

T 

° n^CS ,0) A(T(S 0)) 
m m' 


(26) 


(27) 


is the Initial radiative cooling time scale of the hottest plasma. 

A two parameter family of solutions for (26) can be obtained by setting 

the two time dependent factors to constants. Three equations result: two 

coupled first order i=qMations which can be Integrated completely to yield 

the time dependence of the temperature and pressure (i.e., 0(t) and (j>(t)), 

and one second order equation which can be integrated once to yield the 

d£ 

differential emission measure (i.e., as a function of 5). The solutions 
to (26) are: 


0(t) = (»1 + nt)”^ (27a) 

-V-.1 


(j)(t) = (1 + rit) 


(27b) 
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5 

ds 


(1 + (0) . (e + 

d» 


(28) 


where y and V are the two arbitrary parameters, and g and n are given by: 


b(v) 


3(1 + V) A + 3 

+ ' ■ 

(2£-l)v + 2Jl + 4 2 + 2 


X- V/ 


and 


3(2 + 2) 

n 

((22-l)v + 22 + 4)y T 

o 


(30) 


Equation (28) can be integrated (at least numerically) to yield the temper- 
ature profile; however, we are more interested in expressing the emission 
measure profile as a function of temperature since that is the observed 

quantity. Letting T be the spatial dependence of the temperature, 

s 


1 /2+z 


Tg = ?(s) 


(31) 


yields: 


2 , 2 

n“(0,0) 2+v“'^ . /T 

Q(T,t) . 0(t) . 1 


(32) 


dX (0) 
s ' 

ds 


1 + y 


Equation (34) is the main result of the isobaric model. We wish to 


determine whether a steep dependence of Q on T can be obtained for par- 
ticular values of the parameters y and v. It is evident from (32) that the 
steepest dependence is obtained if 
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Ip I « and g(V) < 0 (33) 

Since in that case, Q ~ T ^ , But, it is shown below that (33) 

8 

implies physically unacceptable solutions. 

Using (22), the velocity can be related to the temperature gradient; 

viz, 

dT 3 

V — ~ — 0' T (T + y) (34) 

ds 5T^y ® ® 

In crder that the velocity be directed opposite to the temperature 
gradient; i.e. , material moves downward as observed, we require, 

dT 

V — < 0 (35) 

ds 

Now, assuming that |y| < T implies we must have y>0 so that (35) is 

s 

consistant with (34) . However, we also require that the temperature and 
pressure, 0(t) and (j)(t), decrease in time, i.e., the plasma is cooling. 

From (27), (30) and that y Is positive, this requirement implies that either 

V > 0 and (2£-l)v + 25,+4 > 0 (36a) 


or 


V < -1 and (2il-l)V + 2)1+4 < 0 (36b) 


But, if either case (36a) or 36b) is valid, g(v) must be positive, equation 
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(29). Therefore, a physically acceptable solution with |)i| < t is 
not possible unless g > 0, which Is contrary to what is required, (33), 
for a steep profile. 

From the discussion above it is evident that the only solutions of 

Interest are those with (ij| > In that case, Q « T over almost 

s s 

all the range of T . To illustrate this point the dependence of Q on T , 
s s 

i.e. , T + y)“Sj is plotted in Figure 2 as a function of T for 

the particular parameters: Ji = 1.5, y = -2 with v •= -2, and y = 1 with 

V = 1 near the singular point T = -y. We note that just as in the static 

s 

radiative model, 6 = 2.5 except near the temperature maximum at T =1. 

s 


III. DISCUSSION 


The main result of the previous section is that for radiatively 
dominated cooling of corona and transition-region flare plasma the differ- 
ential emission measure Q “ and, hence, 6 = £ + 1. It appears unlikely 

that this dependence can account for the observations which imply 6^3. 
Calculations of the radiative loss coefficient A(T) originally did indicate 
a large value for f. in the range 10^ S T < 10^ K, e.g. , Cox and Tucker's 
(1969) results indicate A»1.8. However, these authors did not include 
the contribution from Fe ions to the radiative losses. More recent calcu- 
lations (e.g., Tucker and Koren 1971, Raymond et al. 1976) have tended to 
yield increasingly smaller values for H due to the strong Fe lines formed 
in the 10® - lO’’ K temperature interval. One of the most recent results 
gives 6 = .5, Rosner et al. (1978) and, thus 5 ~ 1.5 which is no larger 
than that predicted by conductive cooling models. 

Therefore, we conclude that the proposal of Underwood et al. (1978) 
is incorrect; radiative cooling cannot produce the steep emission measure 
profile observed in some flares. On the other hand, the conductive models 
are also inadequate. There are several possible explanations for the 
discrepancy. One is that the calculations of predicted line Intensities 
and of A(T) are inappropriate for flare plasma; however, this does not 
seem likely. The calculations assume that the plasma is optically thin, 
is in ionization equilibrium and has standard abundances. Although one 
or more of these assumptions may be violated during the primary heating 
phase, they should hold for the long-lived decay phase. 

Another possibility is that flare heating is responsible for the steep 
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profile. However, in that case strong heating would have to continue 
well into the observed decay phase and have a very specific dependence 
on density and temperature in order to reproduce the observations. 

We believe that the most likely explanation (assuming the observations 
are not simply due to instrumental error), is that the models developed 
so far are inappropriate. Various assumptions nave been made; the main 
one being that all plasma velocities are subsonic. This assumption is 
common to both radiative models considered here, and to the conductive 
dominated models, (Antiochos and Sturrock 1976, 1978). Supersonic velocities 
have been observed in cooling flare plasma, but only at low temperatures 
(< 10®'*’ K) , below those for which the emission measure profile exhibits 
a steep dependence (see Underwood et. al. 1978). However, velocities are 
difficult to observe in higher temperature plasma so that large velocities 
may be consistant with the data. We intend to investigate flare decay 
models with supersonic velocities in a subsequent paper, (Antiochos and 
Sturrock 1979), 
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FIGURE CAPTIONS 

1. Graph of log Q versus log (T(s,t)/T^) for the value of Jl » 1.5. The 
solid curve refers to the time t = 0, i.e. the initial emission measure 
profile, and the broken curve to the later time, t » tj ♦ 

2. Graph of the spatial dependence of log Q, i.t. log [T (T + p)"8], 

8 8 

versus log T for two isobaric radiative models with Si = 1.5. The solid 
s 

curve refers to the case with y = 1 and V « 1, and the broken to the 
case with |J = -2 and v = -2. 


log 0 



Figure 1 






